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存在 非 加 热 起 始 段 的 热 平 板 着 火 分 析 
李 尚 胶 卓 建 坤 唐 勇 姚 强 


(清华 大 学 热 科学 与 动力 工程 教育 部 重点 实验 室 ， 北 京 100084) 


摘 要 对 存在 非 加 热 起 始 段 的 热 平 板 着 火 问 题 进行 了 理论 分 析 , 推导 出 具有 局 部 相似 性 的 解 ， 并 采用 基于 高 活化 能 
设 的 匹配 渐进 展开 法 进一步 得 到 了 理论 解 . 非 加 热 的 起 始 段 可 以 降低 气流 通过 加 热 段 时 壁面 附近 的 流速 , 从 而 更 容易 引发 
着 火 . 结果 表明 ， 当 常温 起 始 段 长 度 LX 和 一 般 的 恒温 热 平 板 着 火 距离 xz 比值 xI/Lx = O (e), 此 时 着 火 距离 XLI 和 Lx 
的 平方 根 成 反比 . 通过 数值 计算 发 现 ， 当 XI/Zx < 3， 局 部 相似 解 的 误差 不 超过 5%。 本文 的 分 析 方 法 可 以 进一步 应 用 到 
二 维 棉 形 体 、 三 维 旋转 体 边界 层 等 更 为 复杂 边界 条 件 下 的 着 火 问题 . 
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Thermal-ignition Analysis in a Boundary Layer on a Hot Plate with an 
Unheated Inlet Section 


LI Shang-Peng ZHUO Jian-Kun TANG Yong YAO Qiang 
(Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, 
Tsinghua University, Beijing 100084, China) 
Abstract The problem of thermal-ignition on a boundary-layer at a hot plate with an unheated inlet 
has been considered, and a local similar solution is derived. An explicit expression for the minimum 
ignition distance has been derived by using matched asymptotic analysis in the realistic limit of large 
activation energy. The unheated inlet section can reduce the flow velocity at the boundary-layer 
of the inlet of heating section, and then lead to easily ignited. It is shown that when the ratio 
between the length of unheated inlet section L\ and the ignition distance on a isothermal hot Plate 
XI is large enough, the ignition distance XLI is inversely proportional to the square root of Lx. 
According to numerical results, the difference between the local similar solutions is less than 5% when 
Xi/Lx < 3. Although attention in this paper is confined to the fat-plate problem, it should be noted 
that this analytical method can be applied to more general cases, such as wedge flow ignition problems. 
Key words thermal ignition of hot surface; unheated starting length; local similar; asymptotic 


analysis 
符 号 表 

v 化 学 计量 系数 ; a 在 式 (22) 中 定义 的 量 ; 
WwW 摩尔 质量 ; hi 在 式 (19) 中 定义 的 函数 ; 
也 环境 压力 ; 入 入 = Z, 用 于 坐标 变换 ; 
Ro 通用 气体 常数 ， 天 在 式 (21) 中 定义 的 函数 ; 
B 反应 速率 表达 式 中 常数 之 一 ; hm 在 式 (25) 中 定义 的 函数 ; 
W 平均 摩尔 质量 ; 有 在 式 (25) 中 定义 的 函数 ， 
T 温度 或 者 无 量 纲 温度 ; 712 在 式 (26) 中 定义 的 量 ; 
Ea 总 包 反 应 活化 能 ， WwW 在 式 (26) 中 定义 的 函数 ; 
Y 质量 分 数 ; X 广义 着 火 距离 , 在 式 (32) 中 定义 ; 
p 气流 密度 ; K 在 式 (34) 中 定义 的 量 ; 
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u 了 方向 速度 ; Bb 热 扩 散 参 数 ， B= Tw — Too; 
v y 方向 速度 ; 各 摄 动 分 析 中 的 小 量 , e = 72 /Ta; 
w 化 学 反应 速率 , 在 式 (5) 中 定义 ; 7 二 x, 用 于 坐标 变换 ; 
P 无 量 纲 温度 , 了 = cpT/h!; 6 6= x, 用 于 坐标 变换 ; 
L 非 加 热 段 长 度 ; [3 €= f'n); 

L() 算 子 , 在 式 (4) 中 定义 ; 0 热 边界 层 内 侧 反应 区 的 一 级 摄 动 温度 ; 
A 在 式 (6) 中 定义 的 量 ; 和 热 平 板 着 火 问题 中 的 邓 克 尔 数 ， 
a 反应 速率 表达 式 中 温度 指数 或 者 在 式 (15) 中 定义 的 量 ; 上 标 
h 烩 值 ; 0 热力 学 标准 状态 ; 
cp 比 定 压 热 容 ; 二 
s 在 式 (8) 中 定义 的 量 ， BS 
7 在 式 (8) 中 定义 的 量 ; i 第 i 种 组 分 ; 
no 热 值 ; Oo 无 穷 远 ; 

U 主流 速度 ; Ww 壁面 ; 

H 动力 慕 性 系数 ， F 燃料 ; 

人 流 函 数 , 在 式 (9) 中 定义 ; O 氧化 剂 ; 

£ 无 量 纲 流 函 数 , 在 式 (10) 中 定义 ; P 生成 物 ; 

9 在 式 (15) 中 定义 的 函数 ， I 着 火 状态 ; 

711 在 式 (18) 中 定义 的 量 ， L 非 加 热 起 始 段 ; 
X 广义 距离 坐标 ; 

0 二 各 


带 有 化 学 反应 的 边界 层 流动 问题 ， 需 要 求解 高 
度 非 线性 的 偏 微分 方程 组 . 反应 源 项 的 存在 破坏 了 
一 般 无 反应 边界 层 流动 中 的 相似 性 ， 沿 流 线 方向 各 
个 纵 截 面 上 的 物理 量 不 再 具有 简单 的 对 应 关系 。 为 
了 理解 该 问题 中 各 个 参数 之 间 的 相互 作用 ， 减 轻 数 
值 计算 的 负担 , 发 展 出 了 多 种 理论 分 析 方 法 . 

Marblet 采用 级 数 展开 法 分 析 了 层 流 条 件 下 ， 
高 温情 性 气体 和 低温 可 燃气 体 的 混合 层 流 动 中 可 燃 
气体 的 点 燃 以 及 火焰 面 的 发 展 ， 在 忽略 了 沿 流 线 方 
向 的 扩散 后 , 计算 出 从 开始 混合 到 点 火 的 着 火 距离 ; 
而 对 于 着 火 点 之 后 的 火焰 面 发 展 ， 则 考虑 了 沿 流 线 
方向 的 扩散 , 并 应 用 冯 . 卡门 的 积分 方法 对 边界 层 进 
行 积 分 求解 Toong 占 使 用 迭代 法 对 热 平 板 边界 层 
着 火 问题 进行 求解 ， 可 以 通过 低 精 度 的 解 迭 代 得 到 
较 高 精度 的 解 . Lawfal 采用 基于 高 活化 能 假设 的 匹 
配 渐进 分 析 法 , 证 明了 壁面 附近 温度 梯度 为 0 适合 
作为 着 火 判 据 ， 并 由 此 得 到 了 汪 止 边界 层 着 火 的 解 
析 解 . Law 使 用 该 方法 还 求解 了 热 平 板 着 火 向、 一 
维 固体 表面 瞬 态 着 火 加 等 着 火 问题 . 因此 ,匹配 渐 
进 分 析 法 更 能 体现 高 活化 能 化 学 反应 一 般 集 中 在 狭 
案 、 高 温 空间 区 域 的 物理 事实 , 求解 过 程 相对 简单， 

上 述 炽热 表面 着 火 问 题 的 理论 解 都 是 建立 在 热 
表面 温度 为 均一 恒定 的 假设 之 上 ， 然 而 在 实际 的 工 
程 应 用 中 ， 炽 热 表 面 着 火 中 的 热 表 面 并 不 是 恒定 温 
度 的 ， 可 燃气 流 在 经 过 热 表 面 之 前 ,往往 已 经 有 了 
一 定 的 速度 边界 层 厚度 。 即 使 对 于 较为 简单 的 热 平 


板 着 火 实验 ,也 难以 让 热 边界 层 和 速度 边界 层 从 同 
一 个 位 置 开始 发 展 , 因此 , 研究 人 员 在 平板 加 热 段 的 
前 缘 设 置 了 速度 边界 层 的 抽 吸 口 57, 不 仅 增加 了 实 
验 闭 置 的 复杂 度 , 还 可 能 对 流 场 产生 干扰 . 所 以 , 有 
必要 研究 存在 常温 非 加 热 的 起 始 段 对 炽热 表面 着 火 
过 程 的 影响 。 

本 文 基于 一 般 形 式 的 边界 层 反 应 流 偏 微分 方程 
组 ， 推 导出 了 热 边 界 层 发 展 滞后 于 速度 边界 层 时 的 
局 部 相似 解 ， 进一步 采用 基于 高 活化 能 假设 的 匹配 
渐 近 展开 法 得 到 了 此 时 着 火 距离 的 理论 解 .最 后 利 
用 数值 计算 对 局 部 相似 解 的 准确 性 进行 验证 。 


1 基本 控制 方程 


如 图 1 所 示 ,， 本文 所 研究 的 物理 问题 是 : 层 流 
条 件 下 , 一 定 配 比 的 可 燃 预 混 气 以 恒定 的 速度 , 通 


速度 边界 层 ” 热 边界 层 


图 1 存在 非 加 热 起 始 段 热 平 板 着 火 示 意图 
Fig. 1 Schematic of thermal-ignition by an isothermal flat 
plate with an unheated inlet section 


3 其 李 尚 鹏 等: 存在 非 加 热 起 始 段 的 热 平 板 着 火 分 析 659 


过 无 限 流域 中 的 一 块 平板 形成 边界 层 反 应 流 ， 平 板 
不 可 渗透 且 无 催化 作用 ,平板 存在 与 来 流 气 温 相 同 
的 非 加 热 起 始 段 , 剩余 部 分 为 恒定 温度 的 高 温 段 . 当 
平板 的 加 热 段 温度 足够 高 、 气 流通 过 的 距离 足够 长 ， 
在 平板 某 个 位 置 的 附近 将 会 引发 着 火 , 产生 火焰 . 本 
文采 用 稳 态 的 方法 分 析 着 火 的 临界 状态 ， 并 不 考虑 
着 火 之 后 预 混 火焰 在 边界 层 中 的 传播 . 

为 了 简化 分 析 过 程 , 考虑 一 步 总 包 不 可 逆反 应 ， 
只 有 两 种 反应 物 , 分 别 是 氧化 剂 [O] 和 燃料 下， 只 
有 一 种 产物 [P], 化 学 反应 式 为 : 


vo [O] + ve [F] — vp [P| (1) 


根据 阿 累 尼 乌 斯 定律 ， 各 个 反应 物 的 反应 速率 
正比 于 : 


W WW ZO 十 ZF 
O 下 


Ea VO WUF 
(- 褒 ) YeoY (2) 
其 中 , B 和 a 和 皆 是 常数 ,，R? 是 通用 气体 常数 , Wi 是 
第 i 种 组 分 的 摩尔 质量 ,vi 是 第 i 种 组 分 的 化 学 计 
量 系 数 ,Y; 是 反应 流 中 第 i 种 组 分 的 质量 分 数 ，E。 


是 该 总 包 反应 的 活化 能 ， 
采用 Law 革 中 的 假设 条 件 ， 可 以 推导 出 边界 层 
反应 流 基本 方程 组 为 : 
Opu Opv 
2 =0 
L(u)=0 (3) 
L (YF) = 一 w 
L(T)=w 


其 中 无 量 纲 化 的 温度 全 记 为 了 = cpT/ho, 在 不 
引起 歧义 的 情况 下 仍旧 使 用 原 有 符号 Tp 是 反应 流 
的 密度 , 算 子 L(-) 定义 为 : 


0 O 0 0 
L()= (mm 六 十 和 训 ( 训 ) (4) 


方程 组 (3) 里 能 量 方 程 、 组 分 方程 中 反应 源 项 
的 表达 式 为 : 

w (T, Yo, YF) 

ho Qa—(votrve) 
= AT®- (Yotrs) (=) exp 

Cp 
ba VOWVF 和 

(六 )W ;0 


Wevr pW Se 
VTS \ RI 


Wovo 
WEevF 


Ws 0 
其 中 ，c 是 混合 流 的 比 定 压 热 容 , 已 经 假设 了 各 个 
组 分 的 比 定 压 热 容 均 相 等 且 为 定 值 , p 是 环境 压强 ， 
是 @ 第 i 种 组 分 的 生成 痊 ，ho 相当 于 该 燃料 的 
热 值 . 

热 平 板 前 面 存 在 常温 起 始 段 , 长 度 为 工 , 将 坐标 
系 的 原点 设 在 加 热 段 进口 ,此 时 的 边界 条 件 如 下 : 


h 风 二 有 网 十 ho 一 


T=—L: w=U, v=0, Yr = Yr,w 
OYFE 
V 一 二 三 VU, Tr———— 
Oy (7) 
一 上 <Z<0， y=0: 了 工 =7- 
2 过 0,V 王 T= 171w 


可 以 发 现 动量 方程 的 边界 条 件 仍然 具有 一 致 
性 ， 可 以 沿用 典型 的 Howarth-Dorodnitsyn 坐标 变 
换 (z,y) 一 (s,n) 让 动量 方程 与 能 量 方程 解 厢 ， 其 
中 : 


s=UpoHw (T+L) 


(8) 
VY i (Z 十 卫 ) ee 人 pz ) 


是 流通 数 , 令 


由 此 可 以 自动 满足 可 压缩 的 连续 性 方程 . 
在 新 坐标 系 下 的 无 基 纲 流 函 数 f (s,m) 为 : 


Up 
f(s,7) = A 


Chapmanls| 基于 实际 气体 的 热 物理 性 质 给 出 一 
个 重要 假设 ,可 以 简化 二 阶 偏 微分 方程 组 中 的 二 阶 
扩散 项 : 


(10) 


Ph = pooloo (11) 


由 此 , 方程 组 (3) 中 的 动量 方程 可 以 化 为 典型 
布 劳 休 斯 方程 的 形式 : 
d3 df 
d= 


f (0) = 0,7"(0),f"(0%)=1 (12) 


然而 由 于 温度 方程 的 边界 条 件 并 不 是 一 致 的 ， 
温度 方程 和 组 分 方程 不 能 通过 上 述 变 换 直接 得 到 相 
似 解 . 
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0 区 域 , 通过 待定 函数 法 进行 计算 , 利用 如 下 的 自 变 
基 1, 通过 对 不 同 处 x 的 y 坐标 进行 玉 度 变换 ， 得 


2 局 部 相似 解 

将 无 基 纲 流 函 数 在 壁面 处 展开 ， /() = 
了 (OP (0) + 0(77),， 计算 4 时 只 取 了 的 
第 一 个 非 零 项 , 即 假设 热 边 界 层 中 沿 平板 方向 的 速 
度 v 与 7 线性 相关 (在 不 可 压 流动 中 ， 相当 于 物理 
空间 中 的 速度 “在 y 方向 是 线性 分 布 的 ): 


d qf 
UD) Tho 


u=U (zx) (13) 

该 假设 的 合理 性 源 自 三 个 方面 : 

1) 由 于 非 加 热 起 始 段 工 ， 热 边界 层 厚度 小 于 速 
度 边界 层 ; 

2) 一 般 的 燃烧 反应 具有 较 高 的 活化 能 ， 
中 在 热 边界 层 底层 ; 

3) 由 布 劳 修 斯 方程 可 以 知道 (3)(0) = 
f 史 (0) = 0, 这 使 得 式 (13) 的 泰勒 展开 式 具 有 三 阶 
精度 . 

在 式 (13) 中 带 入 7 的 具体 表达 式 (8), 此 时 有 


= pUf"(0)4 
pu= pUf"( Tr) p(T,Y ) 


反应 集 


Qg(7,Y) X,Yy) 
et 9 
求解 三 阶 常 微分 方程 (12) 可 以 得 到 上 式 
中 ,，f”(0) s 0.46960, 另外 ， 
=U7"(0) /3—— 
Foe Hee (15) 


y 
oo ) dy’ 
g(x,y) ,| p(z,Yy) dy 
结合 方程 组 (3) 中 的 连续 方程 ， 有 : 


y / / 
m=-| EE) = 
0 TX 


Qa 7 Og / CQ 提 / / 
_ < d 
人 re g (zy ) dy 
(16) 
将 式 (14), (16) 带 入 工 算 子 表达 式 (4), 得 到 : 
_ Qg(z,y) 9 EE 小 
L()= SY +(- 二 人 Bd 
a y | 0 0 6 
| re 扩 电 )] 世 -区 (17) 


对 于 本 文 提出 的 热 边界 层 发 展 滞 后 于 速度 边界 
层 的 热 平 板 着 火 问题 , 主要 研究 存在 热 边界 层 的 > > 


到 相似 解 形 式 : 
ap(z) 有 
| DZ ) dy 
其 中 a 为 任意 的 常数 ， 
hi(z)=(L+z) i — L3 


进一步 整理 得 到 : 


op 
L()= -pokoop : 


1 
02 
BR 
az 0 


a2poo Hooh' I VI+L DA 


aF? 0 


其 中 ， 
F(x,Yy) = ee 到 


二 一 
3apcoHpeejiVZ 十 元 Om 


全 | 作 [ E 其 二 J 


38 关 


(18) 


(19) 


(A, 1)， 带 入 (17) 


(21) 


可 以 让 待定 系数 a 取 某 个 定 值 来 进一步 简化 式 


(20): 
s_167"(0)/ vv Vi 
”8V5 (a ) 

此 时 ， 
TN 人 +m 
2836 3 


3 


Om 


将 算 子 (23) 带 入 式 (3) 中 的 温度 方程 ， 


2 pp ,or 
Sr Fn — 12F (x) 


_ 2V367,» (©) 
f° (0)$ pU 


oT 
M1 5 二 


(22) 


使 +3R -1 (LZ+z)i- | 0 (23) 


可 以 得 


(24) 


(25) 


3 期 


能 量 方程 (24) 同时 具有 对 六 的 一 阶 和 二 阶 偏 
导 项 , 可 以 通过 以 下 定义 的 归 一 化 自 变 基 mo 消除 其 
中 的 一 阶 导数 : 

W (m) 


W (m) = 三 eidm,m = 元 (oo (26) 


至 此 ,存在 非 加热 起 始 段 的 热 平 板 着 火 控制 方 
程 组 (3) 可 以 化 成 : 


Fe 站 了 wa 0; 


3 晕 
2736 (W Co) ee ) F(z) 
f" (0)¥ pU 
通过 后 文 的 摄 动 分 析 硒 明 ， 在 着 火 距离 之 前 , 能 
量 方 程 和 组 分 方程 中 的 Max 项 远 小 于 另外 两 项 , 与 
垂直 壁面 方向 的 扩散 相 比 是 三 阶 小 量 ， 因 此 可 进 一 
步 忽 略 简化 , 得 到 的 常 微分 方程 组 形式 如 下 : 


CU 


d2Yp 2336W (co) 及 (z) 

dn2 FY (0)3 pU (28) 
dT 2%36W (co) Eb (2) 

dn2 了 (0)3 pU 

边界 条 件 为 : 
dyYF 

人 d772 ’ (29) 
m=1:Yr= Yr,w, T= To 


而 对 比 不 存在 非 加 热 起 始 段 的 热 平 板 着 火 局 部 
相似 方程 组 由， 


d2 不 27 
= WwW 
dé2 Uf (0 2 
a : (30) 
0 
de? pUf" (0 


可 以 发 现 方程 组 (27) 和 (30) 形式 完全 相同 , 且 
边界 条 件 以 及 上 自 变 基 m2 与 的 5 取 值 范围 相同 , 所 以 
对 于 温度 边界 层 和 速度 边界 层 不 是 同时 开始 发 展 的 
情况 ,着 火 距离 仍然 可 以 得 到 局 部 相似 解 和 由， 通过 
对 比 (27) 和 (30) 右 端 即 可 简单 地 求 得 : 


2336W (00)” Va (ZLD) 村 271 
f" (0)3 pU pUf" (0)’ 


(31) 
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定义 如 下 的 广义 着 火 距离 外， 


ee 2 We ve 
TD) WoW 


W ZO 十 ZF 一 1 po a 一 (Zvo 十 zF) 十 1 
(更 ) 人) 9 


0 
R Cp 


最 终 解 得 : 


XLI = Lx 3 一 (33) 
其 中 3 
i 1 ~ 
36W (00)? f7 (0)3 
1 
~ 1.0405 (34) 


3.30193 . 0.892982 . 0.469603 

所 以 只 要 知道 了 恒温 热 平 板 着 火 距 离 x1, 就 可 
以 由 公式 (33) 直接 计算 出 任意 广义 非 加 热 起 始 段 长 
度 Lx 条 件 下 的 着 火 距离 Xi。 

当 XUZx 久 1 方程 (33) 可 以 化 为 : 


3 
4 kxI 有 4 二 二 
XL,I 名 3Lx ( 乞 ) 一 3 Lx ? (35) 


此 时 的 广义 着 火 距离 Xi 将 与 广义 非 加 热 起 始 
身长 度 二、 的 平方 根 成 反比 , 随 着 非 加 热 起 始 段 长 度 
的 增加 , 着 火 距 离 减 小 的 速度 逐渐 放 绥 ， 
3 渐进 分 析 

至 此 已 将 非 加 热 起 始 段 的 热 平 板 着 火 控制 方程 
组 (3) 简化 成 了 与 无 非 加 热 起 始 段 热 平 板 着 火 控制 
方程 组 外 相同 的 形式 , 因此 , 可 以 使 用 匹配 渐进 分 
析 方 法 进一步 推导 出 理论 解 (也 可 以 认为 式 (33) 的 
唯一 未 知 量 是 xi， 直接 利用 LawI 关于 恒温 平板 着 
火 距离 理论 解 代入 ). 

首先 是 利用 小 量 x = Bm2/e 扩展 边界 层 内 侧 的 
反应 一 扩散 控制 区 域 , 做 坐标 系 转换 (C,&) (7,x)， 
其 中 7 =C=z, 方程 组 (27) 中 的 能 量 方程 化 为 : 

oT 


一 一 2 (w (00) on) Fi (x) Mm ( 


e\’ oT 加 
Om 


Bj) OX 


2736 (W (00) an) FPF (2) /eae\? 
| f" (0)# pU ( ) 2 
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2 
从 上 式 中 看 出 贡 这 一 项 相 比 于 地 项 是 三 阶 
2 
小 量 ,可 以 消去 , 又 因为 反应 主要 存在 于 四 一 0 的 
区 域 , 所 以 进一步 化 为 ， 
0 (37) 
dn2 六 (0)3 pU Bb 
假设 在 反应 内 区 的 组 分 浓度 变化 对 反应 速率 的 
影响 不 大 (主要 考虑 温度 的 影响 ， 了 和 Yo 作为 
定 值 .本 文 的 数值 计算 结果 表明 ， 在 着 火 距离 之 前 
这 个 假设 是 合理 的 ， 在 边界 层 反应 内 区 进行 温度 摄 
动 由 并 带 入 反应 速率 表达 式 , 方程 (37) 可 以 化 为 
d20 人 
dz = = p= 刘 
其 中 0 为 热 边界 层 反应 内 区 的 一 级 摄 动 温 
度 ，47 相当 于 此 物理 过 程 的 邓 克 尔 数 ， 
了 > (Z) 


(38) 


Ar = VY36W (oo 1” (0 A (39) 
其 中 4 为 恒温 热 平 板 着 火 问题 中 的 邓 克 尔 数 ， 
4- 高 Fe 一 (全 
C 
(中 Uf” (0)? 人 


Ta 
ATe (rotrr)tl exp 6 也 ?到 "5 (40) 


当 -< 1, 则 有 


3V6 
= 3 了 


AL = 1 时 正好 对 应 壁面 附近 温度 梯度 为 0, 化 
学 反应 达到 了 自 维持 状态 ， 可 以 由 此 得 到 存在 非 加 
热 起 始 段 时 的 热 平 板 着 火 距离 理论 解 ， 理 论 解 形式 


a (oo) OE (2) A a) 


0.0 02 04 05 0.8 1.0 


(a) T (L,=2x10) 


r~Fhinsyiv 售 人 
IadAIVD | 


站 hh 
ll 
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上 与 式 (33) 相同 ， 只 是 此 时 x 不 需要 求解 常 微 分 
方程 组 ,有 理论 表达 式 四 


(T, — Te) f° (0)? exp [i ) Ea 
(42) 


XI 二 a 一 (vo 十 zF) 十 3vrz 已 
Te ot ye ye 


4 数值 计算 分 析 

本 章 采 用 二 阶 精 度 的 隐 式 有 限 差分 方法 ， 对 仿 
微分 方程 组 (3) 结合 边界 条 件 (7) 进行 数值 求解. 因 
为 边界 层 方程 属于 抛物 型 偏 微分 方程 ， 可 以 只 对 7 
方向 进行 隐 式 求解 ， 而 z 方向 作为 解 的 推进 方向 . 
数值 计算 采用 的 参数 如 下 四 


a=Vo=Vr=1;Ys,% = Y0,%w = 0.2; 
T, = 103K:;T, = 104K; 人 = 300K 


根据 之 前 的 理论 分 析 ， 壁面 附近 温度 梯度 为 0 
的 点 ( 和 =0) 可 以 作为 着 火 点 .着 火 距离 为 


平板 加 热 段 前 缘 与 着 火 点 间 的 距离 . 文献 [4 恒温 条 
件 下 的 热 平 板 着 火 距离 x1 = 1.56 x 104, 图 2 中 给 出 
的 3 个 算 例 对 应 的 无 量 纲 非 加 热 起 始 段 长 度 L、 分 
别 为 2x104 4x104 和 6x104 与 常温 热 平 板 着 火 
距离 比值 分 别 约 为 1.3, 2.6 和 3.8。 

图 2 显示 了 在 不 同 非 加 热 长 度 下 沿 流 线 方向 
平板 边界 层 的 温度 场 和 组 分 场 分 布 ， 从 温度 分 布 图 
2(a), 2(c) 和 2(e) 中 可 以 看 出 , 存在 非 加 热 起 始 段 时 
温度 边界 层 小 于 速度 边界 层 ， 即 使 已 经 假设 普 朗 特 
数 为 1, 两 者 也 不 再 相似 . Lv = 2x104 时 对 应 Xu = 
1.33 x 104, Lx = 4 x 104 时 对 应 Xu = 1.11 x 104， 
而 Lx = 6 x 104 时 xrr = 0.96 x 104, 对 于 7 不 
是 远 小 于 x 的 情况 , 增加 非 加 热 起 始 段 长 度 使 着 火 


(b) 五 (=2x109 


甘 日 干 


只 月] 
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(c) T (L,=6x10%) 
图 2 不 同 非 加 热 起 始 段 长 度 对 温度 、 组 分 分 布 的 影响 


Fig. 2 The effect of the length Lx of unheated section on the Profiles of temperature and concentration along streamlines 


距离 减 小 。 这 是 因为 非 加 热 起 始 段 长 度 加 长 会 使 热 
边界 层 内 的 流速 相对 减 小 , 从 拉 格 朗 日 坐标 系 来 看 ， 
流体 微 团 在 移动 固定 距离 时 耗 时 更 长 ， 反 应 累积 放 
热量 更 多 , 这 最 终 导致 着 火 距离 缩短 . 

图 3 中 实心 点 代表 的 是 从 不 同 LX 的 热 平 板 平 
板 着 火 数 值 算 例 中 提取 出 的 着 火 距离 ， 而 实 线 是 利 
用 局 部 相似 解 (33) 计算 得 到 的 . 局 部 相似 解 只 需要 
求解 一 次 LX = 0 时 的 常 微分 方程 组 ， 对 于 任意 的 
Lx 都 可 以 直接 得 到 对 应 的 着 火 距离 xL,1。 发 现 当 
Lx > 5000， 也 相当 于 是 Xi/Lx < 3 的 时 候 ， 局 部 
相似 解 和 数值 解 的 误差 在 5% 以 内 , 两 条 曲线 的 趋势 
一 致 ， 可 见 局 部 相似 解 非常 好 地 描述 了 存在 非 加 热 
起 始 段 下 的 着 火 过程 。 局 部 相似 解 之 所 以 在 XI/Zx 
较 大 的 时 候 失 效 ， 是 因为 推导 过 程 中 假设 了 在 加 热 
段 的 附 面 边 界 层 中 无 基 纲 流 函 数 为 线性 分 布 ， 而 这 
显然 在 非 加 热 起 始 段 较 短 的 情况 下 误差 是 比较 大 的 ， 
此 时 速度 边界 层 还 比较 薄 ， 线 性 假设 会 使 得 边界 层 
内 的 流速 被 高 舍 ， 从 而 也 就 高 舍 了 着 火 距离 。 局 部 
相似 解 在 所 有 的 区 域 都 会 稍微 大 于 数值 解 . 

而 对 于 用 匹配 渐进 分 析 法 得 到 的 理论 解 ， 参考 
Law 的 数值 分 析 结 果 四 ,其 误差 源 自 局 部 相似 假设 
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和 匹配 近似 两 方面 ,根据 上 面 的 数值 分 析 ， 可 以 知 
道 局 部 相似 假设 带 来 的 误差 远 小 于 匹配 近似 带 来 的 
误差 (50% 左 右 ), 所 以 理论 解 误 差 主要 源 自 匹配 近似 
的 过 程 。 


一 一 Numerical 
1 一 一 Theoretical 
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图 3 不 同 长 度 非 加 热 起 始 段 下 着 火 距离 


Fig. 3 Dimensionless ignition distance XL,I for different 
unheated inlet length Lx 


5 结 论 


本 文 对 存在 常温 非 加 热 起 始 段 的 热 平 板 着 火 问 
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题 展开 了 求解 ， 基 于 热 边 界 层 中 沿 平板 方向 速度 4 
与 9 线性 相关 的 假设 , 对 物理 坐标 采用 适当 的 尺度 
变换 , 可 以 推导 出 该 问题 的 局 部 相似 解 . 关于 速度 v 
与 n 线性 相关 的 假设 合理 性 源 自 三 个 方面 : 1) 热 边 
界 层 滞后 速度 边界 层 ; 2) 高 活化 能 反应 集中 于 高 温 
区 域 ; 3) 壁面 附近 /3) (0) = /79 (0) = 0. 根据 数值 
计算 , 局 部 相似 解 与 数值 解 吻合 良好 . 

利用 匹配 渐进 展开 法 可 以 进一步 将 局 部 相似 解 
推出 理论 解 ， 理 论 解 的 误差 主要 源 自 匹配 近似 的 
过 程 . 

具体 的 理论 推导 过 程 可 以 进一步 推广 到 小 角度 
二 维 棉 形 体 、 小 角度 旋转 体 等 存在 非 加 热 起 始 段 条 
件 下 炽热 表面 的 着 火 问 题 ， 对 任意 温度 分 布 的 热 平 
板 着 火 问 题 具有 借鉴 意义 。 另外， 由 于 存在 非 加热 
起 始 段 条 件 下 热 边界 层 和 速度 边界 层 两 者 间 本 身 就 
不 存在 相似 性 ， 可 以 去 除 质量 、 动 量 和 能 量 扩散 系 
数 都 相等 的 假设 . 
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